Particle Heating by Alfvenic Turbulence in Hot Accretion 

Flows 

Eliot QuataertQ 

Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 

ABSTRACT 

Recent work on Alfvenic turbulence by Goldreich & Sridhar (1995; GS) 
suggests that the energy cascades almost entirely perpendicular to the local 
magnetic field. As a result, the cyclotron resonance is unimportant in dis- 
sipating the turbulent energy. Motivated by the GS cascade, we calculate 
the linear collisionless dissipation of Alfven waves with frequencies much less 
than the proton cyclotron frequency, but with perpendicular wavelengths of 
order the Larmor radius of thermal protons. In plasmas appropriate to hot 
accretion flows (proton temperature ^ electron temperature) the dissipated 
Alfven wave energy primarily heats the protons. For a plasma with (3^5, 
however, where f3 is the ratio of the gas pressure to the magnetic pressure, 
the MHD assumptions utilized in the GS analysis break down before most 
of the energy in Alfven waves is dissipated; how the cascade then proceeds 
is unclear. 

Hot accretion flows, such as advection dominated accretion flows (ADAFs), 
are expected to contain significant levels of MHD turbulence. This work sug- 
gests that, for /3 ^ 5, the Alfvenic component of such turbulence primarily 
heats the protons. Significant proton heating is required for the viability of 
ADAF models. We contrast our results on particle heating in ADAFs with 
recent work by Bisnovatyi-Kogan & Lovelace (1997). 

Subject headings: accretion - hydromagnetics - plasmas - turbulence 



1. Introduction 

In astrophysical accretion flows with significant angular momentum, the accreting 
gas is believed to form either an optically thick, geometrically thin, disk (Shakura & 
Sunyaev 1973; see Frank, King, & Raine 1992 for a review) or an optically thin, geomet- 
rically thick, quasi-spherical flow (Shapiro, Lightman, & Eardley 1976; Ichimaru 1977; 
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Rees et al. 1982; Narayan & Yi 1994, 1995a, 1995b; Abramowicz et al. 1995). In thin 
accretion disks, the gas cools so efficiently that all of the viscously generated energy is 
radiated locally. The resulting low temperature implies that the mean free path due to 
Coulomb collisions is a small fraction of the size of the disk. Consequently, a purely 
fluid description of the accreting gas is reasonable. By contrast, the high temperatures 
in hot accretion flows entail that the (field-free) mean free path is often comparable to 
the size of the accretion flow. Collective plasma effects are thus hkely to be significant 
(Rees et al. 1982). 

In hot accretion flows, the plasma is usually assumed to be two temperature, with 
the ions signiflcantly hotter than the electrons (Shapiro, Lightman, & Eardley 1976). 
Such a temperature difference is only possible if the ions and electrons are thermally 
decoupled. For low accretion rates (m ^ a^, where m is the accretion rate in Eddington 
units and a is the Shakura-Sunyaev viscosity parameter) Coulomb collisions are too 
inefficient to force a one temperature plasma (Rees et al. 1982). One of the outstanding 
plasma physics problems relevant to accretion theory is whether there are collective ef- 
fects which transfer energy from the ions to the electrons on a timescale short compared 
to the inflow time of the gas, thus invalidating the two temperature assumption (e.g., 
Phinney 1981). The only such mechanism that wc arc aware of in the astrophysical lit- 
erature, due to Bcgclman & Chiueh (1988), is probably not efficient enough to eliminate 
the two temperature nature of the flow (Narayan & Yi 1995b). In this paper we thus 
assume that the two-temperature formalism is valid and that the only thermal coupling 
between electrons and ions is due to Coulomb collisions. 

Recently, there has been substantial work (Narayan & Yi 1994, 1995a, 1995b; 
Abramowicz et al. 1995; Chen et al. 1995; Nakamura et al. 1997) on a class of hot opti- 
cally thin accretion solutions flrst discovered by Ichimaru (1977), advection-dominated 
accretion flows (ADAFs). The deflning characteristic of an ADAF is that it is under- 
luminous for its accretion rate (i.e., L <C Mc^). This arises because "standard" ADAF 
models assume that, by virtue of their larger mass, most of the viscously generated 
energy heats the ions. Since only a small fraction of this energy is transferred to the 
electrons via Coulomb collisions, the total energy radiated (almost all by the electrons) 
is much less than the total energy generated by viscosity (Ichimaru 1977; Rees et al. 
1982). The remaining viscously generated energy is stored as thermal energy of the ions 
and is advected onto the central object. 

The assumption that viscosity heats the ions is crucial for the relevance of ADAF 
models. This enables ADAF solutions to exist so long as the electrons and ions are 
thermally decoupled (m ^ a^). If viscosity were to predominantly heat the electrons, 
(optically thin) ADAF solutions would only exist when the electron cooling time is longer 
than the inflow time of the gas, which occurs provided rh ^ 10~^q;^ for synchrotron 
cooling (Mahadevan & Quataert 1997). Thus, in the context of ADAFs, the issue of 
which particles receive the viscous energy acquires particular importance. In this paper. 
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we present a preliminary investigation of this question. In §6.5 we discuss some related 
work by Bisnovatyi-Kogan & Lovelace (1997). We note that an investigation similar 
to ours, but with somewhat different conclusions, was carried out independently by 
Gruzinov (1997). 

We assume that the energy generated by viscosity is initially converted into large 
scale MHD waves. For simplicity, we focus primarily, but not exclusively, on Alfven 
waves. The large scale waves cascade to smaller wavelengths until they are dissipated. 
The relative heating of ions and electrons is determined by which particle species is 
primarily responsible for the dissipation of the waves. Wc focus on linear collisionless 
dissipation mechanisms, as these should be of principle importance in hot accretion flows 
(§6). 

Collisionless dissipation of MHD waves and the back reaction of the dissipated 
energy on the electron and proton distribution functions has been considered extensively 
as a mechanism for accelerating particles in solar flares (e.g., Melrose 1994; Miller & 
Roberts 1995; Miller, LaRosa & Moore 1996). These ideas have also recently been 
applied to particle acceleration in accretion disk corona (Dermer, Miller & Li 1996; Li, 
Kusunose, & Liang 1996; Li & Miller 1997). Following the seminal work of Kraichnan 
(1965), most of these calculations assume that the turbulent cascade is isotropic, that 
is, that the turbulent energy density at any scale depends only on the magnitude of the 
wavevector, and is independent of its direction with respect to the mean magnetic field. 
Recent work on incompressible MHD turbulence by Goldreich & Sridhar (1995; 1997) 
and Sridhar & Goldreich (1994), collectively referred to hereafter as GS, suggests that 
this assumption is inapplicable for Alfvenic turbulence. This is discussed in more detail 
in the next section (§2). 

We use the GS turbulent cascade to motivate the parameter regime (in wavevector 
space) in which we investigate collisionless dissipation of Alfven waves; we also briefiy 
discuss the collisionless dissipation of other MHD modes (§3). In §4 we qualitatively 
discuss the response of the electron and proton distribution functions to the dissipation 
of waves in a nearly perpendicular Alfvenic cascade (we do not, however, solve the quasi- 
linear equations) . Throughout we concentrate on plasmas for which the ion temperature 
is greater than (often much greater than) the electron temperature and the gas pressure 
is comparable to, or greater than, the magnetic pressure, a regime rarely explored in 
calculations of the dissipation of MHD waves. In §5 we apply our calculations to the 
GS cascade and address two difficult, but important, questions: (1) to what extent is 
the linear analysis of §3 applicable to the strong Alfvenic cascade developed by GS and 
(2) is the dissipation found in §3 strong enough to dissipate the turbulent energy before 
the MHD assumptions used in the GS analysis are invalid. In §6 we apply our results to 
hot, two temperature, accretion flows, in particular ADAFs, and in §7 we summarize our 
results. We have attempted to make §6, which contains our astrophysical applications, 
comprehensible without a detailed understanding of the plasma physics calculations in 
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§3-§5. Towards this end, Appendix A contains definitions of a number of quantities used 
repeatedly in this paper. 



2. Alfvenic Turbulence 

Energy injected into a fluid/plasma on large spatial scales, if it is unable to dissipate, 
builds up to nonlinear amplitudes and cascades to smaller wavelength (larger wavevector) 
perturbations; this continues until dissipation becomes important and the turbulent 
energy is converted into thermal energy. To investigate the range of wavelengths where 
dissipation occurs (the "dissipation range" of the turbulence), two characteristics of the 
nonlinear cascade are particularly important. The first is the cascade time, i.e., the time 
for nonlinear effects to transfer energy from a wavevector ~ k to a wavevector ~ 2k. 
This determines how rapid the dissipation must be to halt the cascade. The second is 
the path of the cascade in wavevector space. Does it depend only on |k| = k or also on 
the direction? If the dissipation is a function of k (and not just k), as it is for MHD 
modes, this distinction is crucial. 

If there is no preferred direction in the fluid, the turbulence is isotropic, i.e., just 
a function of k. For MHD turbulence, however, the local magnetic field picks out a 
direction and so isotropy is not guaranteed. In fact, numerical simulations have long 
shown that incompressible MHD turbulence is anisotropic, with the energy cascading 
primarily perpendicular to the mean magnetic field (e.g., Shebalin et al. 1983). Incom- 
pressible MHD turbulence corresponds roughly to cascading Alfven waves, since both 
the fast and slow MHD modes are compressive. Recent work on Alfvenic turbulence has 
clarified the nature of this perpendicular cascade (GS; Montgomery & Matthaeus 1995; 
Ng & Bhattacharjee 1996). 

Linear Alfven waves satisfy the dispersion relation uj = VAkz, where u is the mode 
frequency, va is the Alfven speed and kz is the component of the wavevector along 
the mean magnetic field (taken to be along the z-direction) . GS argue that Alfvenic 
turbulence naturally evolves into a "critically balanced" state in which the cascade time 
at a scale k is comparable to the linear wave period at that scale. Furthermore, the 
parallel and perpendicular sizes of a wave at any scale are correlated, with kz ~ k']^^L~^^^, 
where L is the outer scale of the turbulence (the scale on which energy is injected). The 
fluctuating magnetic field strength on any scale is given by Bk ~ Bout{k±L)~^/^ , where 
Bout is the excitation amplitude on the outer scale.Q 



^The GS analysis breaks down if the fluxes of turbulent energy parallel and anti-parallel to the mean 
magnetie field are not equal. This precludes application of their theory to the solar wind, but this 
assumption should be applicable in hot accretion flows. 
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In the next two sections we calculate properties of the linear dissipation of nearly 
perpendicular Alfven waves {k± ^ k^). This calculation is motivated by, but does not 
explicitly utilize, the Alfvenic cascade of GS. Because Alfven "waves" in the GS cascade 
only live for ~ a mode period before nonlinear effects transfer their energy to smaller 
spatial scales, it is somewhat misleading to speak of them as waves (the turbulence is 
strong rather than weak in the turbulence theory sense). We will generally treat this 
notion as unproblematic, but in §5, having discussed the results of the linear analysis, 
we will return to the issue of its applicability to the GS cascade. 



3. Linear Collisionless Dissipation of MHD Waves in Two Temperature 

Plasmas 

3.1. Qualitative Considerations 

Waves in a magnetized plasma can in general have electric and magnetic fields both 
perpendicular and parallel to the mean magnetic field. These fields can strongly effect 
the motion of particles through resonant interactions. This occurs when the frequency 
of the wave, in the frame moving with the particle along the field line, is an integer 
multiple of the particle's cyclotron frequency, 

uj — kzVz = nfl, (1) 

where Vz is the particle's velocity along the magnetic field and Q is the relativistic 
cyclotron frequency (e.g. Melrose 1980). When this condition is satisfied, the particle 
and wave are in phase and the wave can efficiently accelerate the particle. In a collision 
dominated plasma, however, such phase coherence is impossible to maintain. 

In the MHD limit, -C f2p and kpp <^ 1, where pp is the Larmor radius of protons 
with the thermal speed and Qp is the proton cyclotron frequency; in this limit the Alfven 
wave dispersion relation is = kzVA and the resonance condition becomes 

va — Vz = nvA^l^jJ- (2) 

For uo <^ Qp, n ^ resonances in equation (^ can only be satisfied by particles with 
Vz ~ nvA^/oj ^ va- Since the thermal speeds of particles are typically of order the 
Alfven speed, there are a negligible number of such particles; n 7^ resonances are 
consequently unimportant. 

For n = 0, resonance occurs when the wave's phase speed along the field line, 
v\\ = uj/kz, equals Vz- Particles with Vz ^ v\\ are accelerated by the wave, while those 
with Vz ^ v\\ are decelerated. Thus, the wave is damped provided that the slope of 
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the particle distribution function at = f y is negative - as it is for a Maxwellian. A 
necessary (but not sufficient) condition for strong damping is that f y be comparable to 
the thermal speed of the particles, so that there are a large number of resonant particles. 

The n = resonance actually corresponds to two physically distinct wave-particle 
interactions. In Landau damping (LD), particle acceleration is due to the longitudinal 
electric field perturbation of a wave (i.e., the usual electrostatic force, E^). In transit- 
time damping (TTD), the magnetic analogue of LD, the interaction is between the 
particle's effective magnetic moment (yU = mv\/2B) and the wave's longitudinal mag- 
netic field perturbation, (Stix 1992). TTD is thus analogous to Fermi acceleration 
(for discussions of the relationship between the two, see Achterberg 1981; Miller 1991). 

In the MHD limit, which corresponds to small wavevectors, the Alfven wave has 
both = and -B^ = and so is undamped by linear coUisionless effects. For larger 
wavevectors, the MHD approximations are less applicable and kinetic theory corrections 
to Ez and become important, leading to finite dissipation of the Alfven wave; the 
n 7^ resonances may become important, if the wavevector has a significant component 
parallel to the background magnetic field. For the perpendicular cascade of Alfven waves 
due to GS, however, when k^pp ~ 1, ~ Vtp{pp/ LY^^ <^ Vlp and so n 7^ resonances can 
be satisfied only by particles with ^ va, of which there are a negligible number. Since 
we expect significant E^ and/or B^, and thus significant dissipation, when kj_pp ~ 1, 
this implies that, even in the dissipation range, only n = resonances are important. 

We now investigate in detail the coUisionless dissipation of Alfven waves, focusing on 
k±pp ~ 1 and k^pp <^ 1. This problem has been considered by a number of authors (e.g., 
Akhiezer et al. 1975; Hasegawa & Chen 1976; Stefant 1976), but not in the parameter 
regime of interest to us. 



3.2. Detailed Calculations: Methods 



Consider a coUisionless hydrogen plasma which is homogeneous, fully ionized, and 
threaded by a mean magnetic field, B = Bqz. We assume that, in the unperturbed 
state, each particle species in the plasma (electrons and protons) has an isotropic, non- 
relativistic, thermal distribution function with no bulk (average) velocities.^] Small am- 
plitude perturbations to the equilibrium state of the plasma satisfy the following dis- 
persion relation, which is obtained by linearizing and Fourier transforming (in time and 
space) Maxwell's equations (Stix 1992; Chapter 1) 

k X (k X E) + — e ■ E = 0. (3) 



^By isotropic we mean that the temperature is the same perpendicular and parallel to z. 
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E is the electric field perturbation and the dielectric tensor, e, is given by e^j = 6ij + 
J2sXij^ where the sum is over the susceptibility tensor (xtj) of each particle species in 
the plasma; for our case, s takes on two values, e and p, for electrons and protons, 
respectively. The susceptibility tensor is calculated by combining the linearized and 
Fourier transformed versions of Maxwell's equations and the collisionless Boltzmann 
equation, and is given by, taking k in the x-z plane, (Stix 1992; Chapter 10) 



, _ mpVtp exp(-A,) ^ , 

Xxx 2 ~ \ J-n'^o\'in) 



/(.xy Xyx 



^pVtp exp(— As) ^ .n^ 



\qs\ Vj^ UJrjAs n=-oo 



Xyz ^zy 



r 2 ~ 2^ --'„;[-'-+ 'Sn^oUr; 



UJT] 



where Vts = {2kBTs/msY^'^, Qs, nis, and T, are the thermal speed, charge, mass, and tem- 
perature of the particles, respectively, va = Bq/ {A-k pY^"^ is the Alfven speed and /„ is the 
modified Bessel function with argument A^ = k'\y^g/2Vll = 0.5k\pl, where ps = Vts/^s 
is the Larmor radius of particles with the thermal velocity and Qg = QsBo/msC is the 
non-relativistic cyclotron frequency (taken to be a signed quantity). Prime denotes dif- 
ferentiation with respect to Ag. Note that we are primarily interested in two temperature 
plasmas so that Tp and Tg are, in general, not equal. 

In deference to the problem of interest, we measure mode wavelengths and fre- 
quencies in terms of the proton Larmor radius and cyclotron frequency, using rj = k^Pp 
and \p = 0.5k'^Pp for the parallel and perpendicular components of the wavevector and 
oj = uj/flp for the mode frequency. We emphasize, however, that equations (0) and (|) 
are valid regardless of the magnitude of rj, Ap, or uj, and can be used to investigate the 
properties of any plasma waves (subject to the validity of the assumptions stated at the 
beginning of this subsection). 

In equation (^, Zg, the plasma dispersion function with argument 

e„ = (5) 
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is given by, taking fc^ > 0, (Stix 1992; Chapter 8) 



where the contour, F, is such that the pole at z = ^ hes above the contour of integration 
in the complex z plane. Mathematically, coUisionless dissipation arises from the con- 
tribution of this pole to the susceptibility tensor, which is oc exp (— by the residue 
theorem. Thus, a necessary (but, again, not sufficient) condition for strong damping is 
^„ ^ 1, which is the thermal average of the single particle resonance condition, equation 

Solving the dispersion relation yields, for a given real k, a complex mode frequency, 
whose imaginary part, 7, represents the growth or dissipation of the wave. It does 
not, however, directly reveal the relative energy absorbed by the protons and electrons 
as the wave is damped. Strictly speaking, this can only be obtained from a nonlinear 
theory since energy considerations are necessarily second order in the amplitude. The 
Ohmic heating law, however, gives that the rate of change of energy of particle species 
s is oc js ■ E, where js is the current. The complication is that the Ohmic heating law 
includes both the increasing thermal energy of the resonant particles (responsible for 
the damping) and the decreasing oscillation energy of the non-resonant particles (e.g., 
Barnes 1968a). When considering heating of the plasma by wave dissipation, one wishes 
to calculate only the former, not the latter. In the weak damping limit (7/co' <^ 1), this 
can be obtained from the Ohmic heating law by relating js to E using the susceptibility 
X evaluated at real frequencies. This yields (Barnes 1968b; Stix 1992) 

_ E* ■ Xs Ilm(<..)=0 • E 

where Ps is the energy absorbed in a mode period, per unit wave energy, by particle 
species s and x1 = iXs — xl)/2^ is the antihermitian part of the susceptibility tensor. 
Physically, Xs is evaluated at real frequencies since this entails that its only contribution 
is from the imaginary part of Zo{C,n)', in turn, for Im(ci;) = 0, the only contribution to 
Im(Zo) is from the poles at 2; = As discussed below equation (1), the contributions 
from the poles in Zq correspond to the thermal average of the single particle resonance 
condition (eq. |]1|). Setting lm{u) = in equation (0) therefore isolates the contri- 
bution from the resonant interactions, which is precisely what one wishes to do when 
determining particle heating. In equation (|^), W, the wave energy, is given by 

w^^mc^w.^i..,).^ (8) 

where e/i = (e + e^)/2 is the hermitian part of the dielectric tensor and B is the wave's 
magnetic field perturbation. We note that, for 'y/u <^ 1, energy conservation implies 
that Pp + Pg = 27T, where T = 27r/Re(co') is the mode period. For the dissipation of 
Alfven waves of interest to us, this is well satisfied even if 7T ~ 1; we therefore always 
use Pp and Pg as estimates of the proton and electron heating rates, respectively. 
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3.3. Alfven Waves with large k± 

Because equation (^) has an infinite number of roots (mostly strongly damped waves 
with no fiuid counterparts), some care must be taken in its solution. Our technique is 
to first solve the dispersion relation in a simple limit (e.g., the MHD limit and plasma 
parameters such that the wave is weakly damped); this ensures that we know which 
mode we are investigating. We then incrementally change the wavevector and/or the 
plasma parameters and follow the properties of the solution. We emphasize that we 
have used the exact form of the susceptibility tensor, with no approximations (save for 
terminating the sum over Bessel functions at some appropriately large number). 

In Figure 1 we show several properties of the Alfven wave as a function of Tp/Te 
and the parameter Ap = 0.5kj_Pp, which measures the perpendicular wavelength of the 
wave; the results correspond to a /3 = 1 plasma, where /5, the ratio of the gas pressure 
to the magnetic pressure, is given hj j3 = SnnkBiTp + T^,)/ = u|,(l + Te/Tp)/v\. In 
the limit of kzPp <^ 1 and Vts,VA <^ c, which we take here, the Alfven wave properties 
given in Figure 1 are independent of the exact values of k^Pp and Vts,VA, and depend 
only on Ap, Tp/T^, and p. 

Figure la shows the parallel phase speed of the Alfven wave in units of the Alfven 
speed {v\\/va = Re(co')/A;2f^). The two curves correspond to Tp = and Tp = lO^Tg. 
For Tp ^ lOTg, f II is nearly identical to the Tp = lO^Tg result shown in the figure. In the 
MHD limit (Ap ^ 1), we have v\\ ~ va, the usual Alfven wave dispersion relation, while 
in the Xp^ 1 limit ~ vaV^p — VAk±Pp, i.e., the wave frequency depends strongly on 
the perpendicular wave number. This result is well-known from analytic treatments in 
the (3 <^ 1 limit (e.g., Hasegawa & Chen 1976). 

The remaining panels in Figure 1 specify properties of the dissipation of the Alfven 
wave by collisionless effects. Figure lb gives the dissipation rate of the mode in units 
of the mode period, 7T. Figure Ic gives the dimensionless proton, Pp, and electron, Pg, 
heating rates, using equation (^. For Tp ^ lOTg, the dimensionless proton heating rate, 
Pp, is nearly identical to the Tp = lO^Tg result shown in the figure. We note that 7T 
(Fig. lb) can be derived from Pg (Fig. Ic) using Pp + P^ = 27T. Finally, in Figure Id 
we explicitly show the relative heating of the protons and electrons, Pp/P^. 

As we are also interested in plasmas with (3^1, Figure 2 shows 7T (Fig. 2a) and 
Pp/Pe (Fig. 2b) for several (3 with Tp/T^ = 100. For /5 ^ 1, the behavior of the dissipation 
of the Alfven wave with varying Tp/T^ is very similar to that shown in Figure 1 (the 
(3 = 1 case): 7 is relatively independent of Tp/T^ for Ap ^ 1, but the relative heating of 
protons and electrons (Pp/Pe) increases with increasing Tp/Te, as in Figure Id. This is 
explained analytically below. 

For our purposes, the relevant parameter regime in wavevector space is Ap ^ 1 . This 
is because, as is discussed in more detail in §5.2, the turbulent dynamics of Alfven waves 
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due to GS is only valid in this limit. For Ap ^ 1, the key qualitative results are that 
the Alfven wave can be strongly damped by coUisionless effects, particularly if /5 ^ 1. 
Furthermore, provided that Tp ^ Tg, the dissipated wave energy primarily heats the 
protons. We note that, for Ap ^ 1, I'll ~ va\J\) ~ Vtp\f\> and so with increasing Ap 
there are progressively fewer protons available to resonate with the wave. This is why 
the proton contribution to the damping, Pp, falls off exponentially for Ap ^ 1 (Fig. Ic). 

The dissipation properties of the Alfven wave shown in Figures 1 and 2 can be 
understood in terms of primarily TTD by the protons; the generally smaller electron 
contribution is due to both LD and TTD for /3 ~ 1, while TTD dominates for larger 

This can be seen by looking at the relative contributions of TTD and LD in the 
expression for the particle heating, P, (eq. 0). The TTD contribution is oc Xyyl-^'yP 
(which is oc Xyyl^zl"^ by Faraday's Law) while the LD term is oc xLl-^^P- Figure 3 
we show, for Ap = 0.1, the ratio of the TTD and LD contributions to the proton (Fig. 
3a) and electron (Fig. 3b) heating rates as a function of Tp/Tg for several /3; below 
we discuss the physical/analytical origin of these results. We emphasize that, even for 
Ap ~ 1, u; <^ fip (since u ~ k^VA and k^Pp <^ 1 by assumption); the magnetic moment 
of a particle is thus an adiabatic invariant, as is required for TTD. 

From equation (^), it is relatively straightforward to find the dependence of Xs '^^ 
(5 and Tp/Tg.Q For a fixed Tp/Tf,, this yields Xyy/xlz ^ P both protons and electrons, 
while for a fixed P we find that xly,e/Xzz,e ^ (J'e/Tp)^ and that xly,p/xtz,p is independent 
of Tp/Te. Physically, increasing /3 (at fixed Tp/Tg) increases the proton and electron 
thermal speeds with respect to the Alfven speed, thus increasing the magnetic moment 
of the thermal particles. This makes TTD more important, which is why Xyy/xlz 1^- 
Increasing Tp/Tg (at fixed /5), on the other hand, effectively decreases the magnetic 
moment of the electrons, which tends to decrease the contribution to electron heating 
from TTD. This is why Xyy,elx1z,e ^ (Te/Tp)^. The parallel electric field of the Alfven 
wave, however, decreases with increasing Tp/T^. In the /5 -C 1 limit, Ez oc Te/Tp (at 
fixed /?), while Ey is independent of T^/Tp (e.g., Hasegawa & Chen 1976; Melrose 1986, 
p. 178)0; we find that this is also roughly satisfied if /3 ^ 1. These scalings imply that 
XlyjEy\''/xlz,p\Ez\'' ~ /?(Tp/Te)2 aud that xly,e\Ey? /xlz,e\Ez\'' ~ /3, which reproduce 
the numerical calculations reasonably well. 

These considerations show that, to a good degree of accuracy, one can take Pp/ Pe — 
Xly,p/xly,eJ since both the electron and proton heating is primarily due to TTD. This 
greatly simplifies evaluating Pp/ P^ analytically since a detailed expression for the electric 



■^For /3 ~ 1 and Tp ~ the Landau and transit time contributions are comparable for both electrons 
and protons. 

"^For the regime of interest here, one need only keep the n — terms and the leading order Ap terms. 

^This is because Ey arises from keeping kinetic terms which are dropped in the MHD limit while 
is due to both thermal and kinetic corrections. 
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field vector is not needed. Using equation (|^) and the MHD Alfven wave dispersion 
relation, we find that 

where we have taken Vte 3> va and Ap ^ 1. The small deviations from equation (|) in 
Figures Id and 2b (for Ap ^ 1) are due to the contribution of LD to the electron heating. 
In the next section we give a more general version of equation which is valid for any 
wave damped primarily by TTD. 

Finally, we note that the results given here for the dissipation of nearly perpen- 
dicular Alfven waves in a /? ^ 1 plasma differ from those obtained by Stefant (1976), 
who found that the Alfven wave dissipation rate reaches a maximum (7T ~ 0.1) at 
some particular value of (3 and decreases for larger 13. The reason for this discrepancy 
is straightforward. Stefant did not solve the full dispersion relation (eq. 0), but in- 
stead used a simplified dispersion relation that neglected all contributions from the yy 
components of the susceptibility tensor.^ This amounts to considering only LD of the 
Alfven wave. This is valid only in the /9 <^ 1 limit when the magnetic compression of 
the Alfven wave, 5^, is negligible (e.g., Hasegawa & Chen 1976, who use a dispersion 
relation very similar to Stefant's, but explicitly state that it is only valid for small j3). 
Our calculations agree with Stefant's in the /3 -C 1 limit; including the yy susceptibilites, 
which are responsible for TTD, is, however, necessary in the /3 ~ 1 limit. In fact, in the 
limit of /5 3> 1, Foote and Kulsrud (1979) have shown analytically, by expanding the sus- 
ceptibility tensor to leading order in that only TTD contributes to the dissipation 
of the Alfven wave. 



(9) 



3.4. General Relations for Particle Heating by TTD and LD 



Equation (^) for the relative heating of protons and electrons in the dissipation 
of a nearly perpendicular Alfven wave is more general than it might appear. From 
the susceptibility tensor (eq. 0]) and the particle heating rate (eq. [0) one can show 
that, for any wave which is damped solely by TTD, the relative heating of protons and 



electrons in the Ap ^ 1 limit is 



P. 



TTD 



rp . 1/2 



exp 




f II 

Vte 



(10) 



^His dispersion relation also assumed Xp <C 1. 
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A similar analysis for LD shows that, for any wave which is damped solely by LD, the 
relative heating of protons and electrons in the Ap ^ 1 limit is 



P. 



e/ LD 





The latter equalities in equations ( JIDD and ( JIT] ) correspond to either subthermal waves 
(f|l ^ Vte,Vtp) or equal electron and proton thermal speeds {vte — Vtp). The physical 
origin of these relations is as follows. The LD expression can be written as {Pp/Pe)LD — 
{me/mp){Vf^/Vfp). The first term {me/mp) is the relative acceleration of protons and 
electrons for a given force. The second term {Vf^/vfp) is the relative number of particles 
available to resonate with the wave (i.e., the relative slopes of the proton and electron 
distribution function at the wave's phase speed, taking t>|| ^ Vye,Vtp)- Equation (p!OD for 
TTD can be written as {Pp/ Pe)TTD — {^e/^p){vf^/'vfp){fip/ f^l), where /i^ oc Tg is the 
magnetic moment of particles with the thermal velocity. The first two terms in the TTD 
expression are identical to the LD expression, and have the same physical interpretation. 
The last term reflects the fact that in TTD, the wave-particle interaction is a function 
of the particle's magnetic moment. The larger Hs, the stronger the coupling between the 
wave and the particle. In LD the corresponding term is q^/q^ = 1 since the wave-particle 
coupling is the electrostatic force. 



Equation (|TOD shows that, quite generally, TTD is a natural mechanism for prefer- 
entially heating protons in plasmas with Tp ^ Tg. This is because in such plasmas the 
protons have the larger magnetic moment and so couple better to a wave's magnetic field 
perturbation. LD, on the other hand, leads to preferential electron heating in plasmas 



with Tp Tg 



3.5. Collisionless dissipation of the fast and slow MHD modes 



A given excitation at the outer scale will, in general, contain both noncompressive 
(Alfvenic) and compressive (fast and slow mode) components. For completeness, we 
therefore briefly consider the collisionless dissipation of the fast and slow MHD modes 
in plasmas with Tp ^ Tg. Since these modes have either Bz or Ez in the MHD 
limit (and as there is no detailed theory of fast or slow mode turbulence to indicate 
if the cascade is parallel, perpendicular, or isotropic) we consider the dissipation only 
in the MHD limit. For Tp ~ Tg this problem has been considered in detail by Barnes 
(1966; 1967; 1968a; 1968b). Figure 4 shows the dimensionless dissipation rate, 7T, of 
the fast (Fig. 4a) and slow (Fig. 4b) MHD modes as a function of 6, the angle between 
the wavevector and the background magnetic field, for several Tp/T^ for a /3 = 1 plasma 
(we take kpp <C 1 and va, Vts <^ c). The Tp = Tg results are identical to those of Barnes 
(1966). 
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It is well known that the fast mode is damped primarily by TTD (Barnes 1966; 
Miller 1991). This is because, in the MHD limit, the electric field vector of the wave 
is along y, which leads to a strong compressional magnetic field perturbation, B^. The 
relative heating of the protons and electrons is thus given by equation ( jlD]) of the previous 
section: 

(S) .M'^.pL(i_p-)(i4),ecU, (12, 

\Pe/ fast ymeTeJ [ \ TeTUp J \ Tp J J 

where we have taken /3 ~ 1 and have used the MHD dispersion relation uj ~ kv^ (which 
is also reproduced by our kinetic theory calculations), so that v\\ ~ sec(^)fA. In a one 
temperature plasma, electrons are preferentially heated by fast modes with 6 ^ 60°. 
This is because, at these angles, Vte ^ f y ^ Vtp and there are no resonant protons, but 
plenty of resonant electrons (see also eq. [|l^). In a Tp ^ Te plasma, however, there 
is no propagation angle for which the fast mode is strongly damped and the electrons 
are preferentially heated. This is because Vte ~ Vtp and so the above condition on the 
wave's parallel phase speed cannot be obtained. This also follows directly from equation 
(|T2|) by setting vtp ~ Vte, in which case {Pp/ Pe) fast — (m^pTp/meTeY^'^ , independent of 
6. In the MHD regime, fast mode turbulence in plasmas with Tp ^ Tg should therefore 
lead to primarily proton heating. There is, however, a large range of propagation angles 
{9 ^ 60° — 70°; see Fig. 4a) for which the fast mode is essentially undamped (as opposed 
to 6* ^ 88° for a one temperature plasma). These modes would likely cascade out of the 
MHD regime. 

The slow MHD mode is essentially a sound wave modified by the presence of a 
magnetic field; it thus has a large Ez and is Landau damped. For Tp ^ T^ this leads 
to preferential electron heating (§3.4). As Figure 4b indicates, the slow mode is very 
strongly damped by collisionless effects; this is particularly true in a Tp ^ Te plasma, 
since there are more electrons available to resonate with the wave. This likely precludes 
slow mode turbulence cascading out of the MHD regime (again, particularly in a Tp 3> Tg 
plasma) . 



4. The Effect of Wave Dissipation in a Perpendicular Alfvenic Cascade on 
the Electron and Proton Distribution Functions 

Particle acceleration by Alfvenic turbulence has often been discussed utilizing ei- 
ther parallel or isotropic turbulent cascades (e.g., for solar flares. Miller & Roberts 1995; 
Smith & Miller 1995; for accreting black holes, Dermer et. al. 1996; Li et. al. 1996). 
In the case of an isotropic cascade, these works also neglect the n = wave particle 
interactions (which are the focus of this paper). These assumptions allow a dramatic 
simpliflcation of the quasi-linear diffusion equations, which describe the response of the 
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particles to the dissipated wave energy as a diffusion in velocity space (e.g., Melrose 
1980). In contrast to the perpendicular cascade of Alfven waves considered here, dis- 
sipation occurs when uj ~ Vtp. This corresponds to the \n\ = 1 resonance in the single 
particle resonance equation (eq. [jT|). In this case the turbulent energy is transferred 
entirely to the protons, regardless of Tp/T^ (since the damping is due to a resonance 
between the proton's cyclotron motion and the wave's perpendicular electric field). 

Dissipation of Alfvenic turbulence by the cyclotron resonance is attractive because it 
naturally leads to the formation of strong non-thermal features in the proton distribution 
function (and acceleration of protons to relativistic energies), which are suggested by 
observations of both solar flares and accreting black holes. To see qualitatively how 
this occurs, note that an isotropic Alfvenic cascade has waves with frequencies from the 
outer scale frequency up to ~ Qp. From the resonance condition, equation (|I|), we see 
that waves with ~ fip accelerate protons with f z ~ f a ~ Vtp (we take /5 ~ 1 in this 
section). Waves with u ^ Qp can only accelerate particles if there is a big Doppler shift 
in the wave frequency, i.e., particles with Vz ~ Vtp^lp/uj. Thus a spectrum of waves with 
uj/Vtp ranging from <^ 1 up to ~ 1 can naturally accelerate particles from thermal to 
relativistic energies. The above references contain detailed calculations of the evolution 
of the proton distribution function by this process. 

The work of GS suggests, however, that isotropic Alfvenic cascades are unlikely to 
be obtained. As we noted in §3.1, when dissipation in a nearly perpendicular Alfvenic 
cascade occurs, u <^ f2p. This implies that n 7^ resonances (such as the cyclotron 
resonance) are unimportant since there are a negligible number of particles with Vz ~ 
Vtpflp/uj ^ Vtp. Since only n = resonances are relevant, the resonance condition 
simplifies to u = kzVz or, using the MHD Alfven wave dispersion relation, f ^ = f a ~ 
Vtp. We emphasize that this condition, which is independent of the wave frequency 
and wavevector, holds for all waves in a nearly perpendicular Alfvenic cascade. The 
effect of the turbulent energy on the proton distribution function is thus qualitatively 
as follows. All waves dissipate their energy to particles with f 2 = f a ~ Vtp. This 
increases the parallel energy of particles near the peak of the thermal distribution, but 
no suprathermal feature is formed; there are simply no waves which can accelerate 
suprathermal particles. Only the parallel energy of the particles increases since, for 
UJ ^ Qp, a particle's magnetic moment is an adiabatic invariant; its perpendicular 
energy therefore cannot change. The resulting parallel proton distribution function is 
reasonably well approximated as thermal, since it is roughly monotonic and has a well- 
defined mean energy (e.g., Begelman & Chiueh 1988). For Tp ^ Tg, the same holds 
for electrons since then vte ~ va and the particles which are heated also reside near the 
peak of the thermal distribution function.]] 



^For Tp ~ Te, the electron heating may lead to a "bump" in the distribution function at Vz — va 
Vte (in addition to the thermal peak near vte)- The same holds for protons if /? is much different from 
1 since then vtp ^ va- In these cases, there is still no acceleration of particles to relativistic energies, 
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The fact that TTD only heats the parallel energy of the particles is significant since 
the efficiency of TTD depends on the particle's magnetic moment, which is proportional 
to the particle's perpendicular energy. If Ej_ is unchanged as the plasma is heated the 
efficiency of TTD would be substantially reduced. In §6.4 we discuss mechanisms which 
lead to isotropy in the distribution function in hot accretion flows. 



5. Applications to the GS cascade 



5.1. Validity of the linear theory 



For the linear analysis of the previous sections to be applicable, a particle's mo- 
tion must be reasonably well approximated by the guiding center approximation (free 
streaming along Bq and cyclotron motion perpendicular to Bq). The presence of a 
turbulent cascade clearly perturbs a particle's motion, so it is worthwhile examining to 
what extent the linear analysis is valid for the strong turbulent cascade of GS. First we 
consider motion parallel to the mean fleld, which is particularly important since we are 
interested in the n = resonances, which occur between the free streaming motion of 
the particle and the wave's parallel flelds. 

For large wavevectors the Alfven wave has flnite perturbed electric and magnetic 
flelds along the mean magnetic fleld {E^ and B^), which are responsible for the damping 
discussed in the previous sections. If the parallel flelds are too large, however, they 
trap the particles in the potential fluctuations of the wave, violating the assumption 
of free streaming. This occurs within the cascade time provided that UTe^t ^ 1, where 
Te ~ {m/ qk^EzY^"^ is the characteristic oscillation period in the potential well of the 
parallel electric fleld (Stix 1992) and ~ [m/ nklBzY^"^ is its analogue for the parallel 
magnetic fleld. From the numerical calculations, E^ ~ E^X^^'^i^T^./Tp and Ey ~ E^r] 
(taking /? ~ 1; see also Melrose 1986, p. 178); using Faraday's Law and the GS scaling 
for By = Bk from §2, 

,^.SJ^^{^" (13) 



and 



Bz^B,.^Al'\^f' . (14) 

Even though the turbulent energy density (S^) decreases with increasing k (§2), the 
parallel fleld strengths increase with increasing k because the kinetic theory correc- 
tions win out over the decreasing wave energy. Using equations (0) and (H), it is 



although a potentially significant non-thcrnial (non-monotonic) feature in the distribution function may 
develop (depending on the ratio of the turbulent energy dissipated to the thermal energy near Vz '^va)- 
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straightforward to check that, for the protons, cur^ ~ X~^^^{L/ ppY^^{Tp/TeY^'^ ^ 1 and 
ujTi, ~ \~^^^{L/ ppY^^ ^ 1 (for the apphcations discussed in §6, L/pp ^ 10®). Thus, 
the turbulence does not appear to contain sufficiently large parallel fields to trap the 
particles. 

We now consider motion perpendicular to the mean field, which is more complicated. 
For a given wave in the GS cascade the E x B drift velocity is ve = cE x B/i?^ ~ 
VAikj_L)^^/^{—y + Ey/Exit). Using the above scahng for Ey/E^ implies 



Ve ~ va 



X 



(15) 



This leads to a Doppler shift in the frequency of the wave seen by the particle of k ■ ve = 
k±VE,x ~ ^'^l^^ipp/Ly^^^ since k is in the x-z plane (§3.2). In addition, the finite lifetime 
of Alfven waves in the GS cascade introduces shifts in the mode frequency which are 
~ u! (since the cascade time is of order the linear mode period). 

Since the Doppler shifts due to the GS cascade are ~ w <^ Qp, the magnetic moment 
of the particles is still an adiabatic invariant. This is important for the applicability of 
TTD. The single particle resonance condition (eq. [|l[]), however, shows up as a delta 
function in the linear theory. It therefore cannot be rigorously applicable to the GS 
cascade since there are additional frequency shifts ~ u not accounted for in the linear 
theory. This leads to a broadening of the resonance condition, but only by a factor of 
order unity. The usual result of resonance broadening is that it makes more particles 
available to resonate with the wave. For our problem, this has little effect since in linear 
theory there are already a significant number of resonant particles - the parallel phase 
speed of the Alfven wave is comparable to the proton and electron thermal velocities. 
Furthermore, while the turbulence is strong in a turbulence theory sense, the wave 
amplitudes are not so large as to broaden the resonance condition by many harmonics 
of the wave frequency, which would significantly reduce the particle heating (Begelman 
& Chiueh 1988). These considerations suggest to us that the linear analysis should be 
a good first approximation for the GS cascade. This does not, of course, preclude that 
other dissipation mechanisms neglected in the linear analysis could be important. 



5.2. Resolution of the GS Cascade 



In the GS cascade, the cascade time is of order the linear Alfven period; conversion 
of a significant fraction of the turbulent energy to thermal energy via dissipation thus 
requires a damping time of order the mode period, i.e., 7T ~ 1. Furthermore, to 



This happens when the Doppler shifts are much greater than the mode frequency. 
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consistently apply the turbulent dynamics of GS, this dissipation must occur when the 
MHD approximations made in their calculation are reasonably applicable. This requires 
kpp ^ 1 and u ^ Qp, which simplify to k±pp ^ 1 for the nearly perpendicular cascade 
of GS. 

From the calculations described in §3, it follows that (for (3^1 plasmas) Alfven 
waves with k±pp ^ 1 and u <^ Qp have damping rates satisfying 7T ^ 1 only if /5 ^ 5. In 
large /? plasmas the proton magnetic moment couples better to the wave's magnetic field 
perturbation, leading to more efficient TTD. If /3 ^ 5, most 50%) of the turbulent 
energy in Alfven waves is dissipated in the GS cascade. The relative heating of protons 
and electrons by this dissipated energy is given approximately by equation (|^). 

What happens, however, for the astrophysically important case of a /5 ~ 1 plasma? 
Proton heating still dominates over electron heating, but Figure 1 gives 7T ~ 0.1 for 
an Alfven wave with k±pp ~ 1, so that only a small fraction of the Alfvenic energy is 
dissipated. TTD is not strong enough to damp the waves before they cascade out of 
the MHD regime. We can envision three possibilities for how the energy is ultimately 
dissipated, but are unable to ascertain which is realized. 

1. Other dissipation mechanisms are important and dissipate the wave energy before 
k± k, pp^. 

2. The cascade continues past k±_ ~ but how it does so (i.e., along what track in 
k space) and in which modes the energy resides is unknown; how the turbulent energy 
is ultimately dissipated is thus unknown. The structure of the GS cascade is crucially 
dependent on the polarization (i.e., the velocity eigenfunction) and dispersion relation of 
the Alfven wave. Both of these properties are entirely different for k^pp ^ 1 (the kinetic 
limit) than they are in the GS regime {k±pp ^ 1). For example, the wave frequency is 
oc k± in the kinetic limit while it is independent of k± in the GS regime. In the GS 
regime, the electrons and protons have the same velocity structure (the E x B drift 
is independent of q and m), while in the kinetic limit, the electrons move substantially 
faster than the protons;^ a single fluid analysis is thus no longer applicable. It is therefore 
unclear whether the turbulent energy will stay in the "Alfven" wave in the kinetic limit. 
Furthermore, GS argue that, in the k±pp ^ 1 limit, the Alfven wave in the nearly 
perpendicular cascade is poorly coupled to the slow and fast modes; this conclusion also 
rests on the polarization and dispersion relations of the modes and is thus inapplicable 
in the kinetic limit. For k±pp ^ 1, the Alfven wave therefore may (or may not) efficiently 
couple to fast or slow (or other) modes. Understanding this possibility is clearly a rather 
complicated problem in plasma turbulence. 

3. As a result of poor coupling to other waves, the energy is dissipated at k±pp ~ 1. As 



^This is because the wave's electric field averages out over the proton Larmor orbit, but not the 
electron's; the protons thus see an effectively smaller electric field and move more slowly as a result. 
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in the previous possibility, the energy "wants" to cascade to higher but being poorly 
coupled to any other modes, it can't; the cascade time and mode energy thus increase 
(maintaining a constant energy flux) until the waves are dissipated. 



6. Applications to Hot, Two Temperature, Accretion Flows 

6.1. General Considerations 

Angular momentum transport in thin accretion disks is now believed to arise from a 
magneto-rotational MHD instability discussed extensively by Balbus and Hawley (1991, 
hereafter BH; for a recent review, see Balbus & Hawley 1997). This instability has been 
considered primarily in the MHD limit, and its applicability to nearly collisionless sys- 
tems is perhaps unclear. In Appendix B we argue (but do not prove) that the instability 
should proceed in collisionless systems, provided that the particle distribution functions 
are close to thermal. The BH instability, when it reaches nonlinear amplitudes, is large 
scale [k ~ H^^ ~ R^^, where k is the wavevector of the instability, H is the disk scale 
height, R is the local radius in the accretion flow, and the latter equality is for quasi- 
spherical accretion flows) and so naturally couples to long wavelength waves, generating 
MHD turbulence; this is seen in numerous numerical simulations (e.g.. Stone et al. 1996). 
This is the basic reason for supposing that most of the gravitational potential energy 
released by viscosity resides in MHD turbulence.[^ In this paper we have made the 
further simplification of focusing primarily on Alfvenic turbulence; this is both because 
Alfvenic turbulence is (comparably) well understood and because some such restriction 
is a necessary first step in attempting to understand the problems considered in this 
paper. MHD turbulence generated by the BH instability may in fact be predominantly 
Alfvenic since the instability is, to linear order, noncompressive. This suggests that the 
excitation of compressive MHD turbulence (such as fast and slow modes) may be less 
important since it is a higher order nonlinear effect than the excitation of noncompressive 
(Alfven) modes.0 

Advection dominated accretion flows (ADAFs) are the only known thermally and 
viscously stable, dynamically consistent models of hot, two temperature, accretion flows 
(Kato et al. 1996). For this reason, we frame the discussion in this section in terms 
of ADAFs. Much of what we say, however, will apply to any hot, two temperature. 



^°More generally, any large scale instability in the plasma will generate significant levels of MHD 
turbulence. 

^^The Alfven wave in a compressible medium is incompressive to linear order, but compressive when 
nonlinear effects are included (Holwegg 1971). 
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accretion flow. For the present purposes, the relevant properties of ADAFs are well 
described by the self similar solution of Narayan & Yi (1994; 1995b), which yields 



^ 2 X 10'' JL_r-^ K, 



-3 



n~ 6.3x10^^ a 'm 'mr cm 

B ^ 10^ ^ ^-1/2^1/2 ^-5/4 q^^^^^^ 

~ lO^^a-l/^ ^ ^-1/2^1/2 ^-5/4 ^g^^ 

^ ~ 6 X 10-V'/'«'/'m-l/2^-l/V-l/^ (16) 

where v/cis the radial velocity in units of the velocity of light, n is the number density of 
electrons/protons, B is the magnetic field strength, determined by assuming a constant 
P in the accretion flowQ, Qp is the proton cyclotron frequency, Pp is the Larmor radius of 
thermal protons, m is the mass of the central object in solar mass units, M = m Mq, rh 
is the accretion rate in Eddington units, M = fn M-^aa (^Edd = 1-39 x lO^^m g s~^), and 
r = R/Rs is the radius in Schwarzschild units {Rs = 2.95 x lO^m cm). In equation ([T6|) , 
the fraction of the viscously dissipated energy that is carried inward by the accreting 
gas is taken to be ~ 1. 

The electron temperature in ADAF models typically saturates at Tg ~ 10^ — 10^° K 
in the inner 10^ — 10^ Schwarzschild radii since the efficient cooling of relativistic electrons 
prevents higher temperatures. Thus the electrons and protons are both marginally 
relativistic and the non-relativistic analysis employed in this paper is a reasonable first 
approximation. Using equation (|16D we can compute the characteristic frequency (z/) and 
(field-free) mean free path {£) for proton-proton (pp) and electron-electron (ee) Coulomb 
collisions (Mahadevan & Quataert 1997): 

z/pp ~ 30q; — - — j m 'm Hz, (17) 

Pee ^ lO^a-^m-^mr^^/^Tg^^/^ Hz, (18) 

^^^^ 

^ ^ lO-W^Vi/^Tl, (20) 
R 

where Tg is the electron temperature in units of 10^ K. Vf.e i^pp) is the rate at which 
the energy and direction of an electron (proton) changes appreciably through Coulomb 



^^The "/3" used in papers on ADAF models, /3adv, is taken to be the ratio of the gas pressure to the 
total pressure and is thus related to the plasma physics /3 used in this paper by /3adv = /3/(/3 + 1). 
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collisions with the same particle species. Electron-proton collisions change the proton 
energy at a rate ~ z/geme/mp. ^ is the mean free path in the absence of a magnetic 
field. Since ^ is often ^ i?, this is indicative of the collisionless nature of the plasma in 
ADAFs. 

The characteristic frequency at the outer scale of the accretion flow, which is roughly 
the frequency at which waves (both Alfven and slow and fast MHD modes) will be 
excited, is Uout — vaR~^ ^ 6 x 10^(/? + l^~'^/'^rn~^r~^^'^ rad s~^. Comparing this with 
I'pp, we see that the protons are effectively collisionless for all perturbations of interest 
(those with frequencies ^ ujout)- By contrast, however, 

Uee/uJout ~ lO^'a-'iP + 1) i/^^Tg"'/' . (21) 

Provided that 

m ^ 10-^a{f3 + l)-^/^T^^\ (22) 

Uee ^ uJout and the electrons must be treated as collisional on the outer scale. This is 
probably not a significant complication for Alfven waves, since collisionless effects are 
unimportant at the outer scale. It must, however, be taken into account in treatments 
of the fast and slow modes, which undergo collisionless dissipation even in the MHD 
limit (§3.5). The fast mode is damped primarily by the protons. The net damping of 
the wave will therefore not be significantly modified by the collisionality of the electrons. 
The slow mode, on the other hand, is a modified sound wave which is strongly Landau 
damped by electrons in a collisionless plasma with Tp ^ Tg. The collisionality of the 
electrons for perturbations with oj ~ Uout will suppress the electron contribution to the 
damping and thus remove this source of electron heating (for accretion rates satisfying 
eq. ||2^). The slow mode will, however, still be strongly damped by other mechanisms 
at the outer scale, e.g.. Landau damping by the protons and wave steepening leading to 
(collisionless) shocks. 

For the GS cascade, strong damping of the Alfven wave occurs when k_[_ ^ ^, at 
which point ~ p~'^^^R~^^^ and lu = uoin ~ Vlp(3~^^'^{pp/ RY^'^ . Using equation (p!6D and 
comparing ujin with z/ge, we find that 

ooju,, ^ Wa"^(3-^'\(3 + ir^'^m^l^m-^l\^'^T^'\ (23) 

Thus, in the dissipation range, collisionless theory can be consistently applied for both 
electrons and protons. This conclusion is strengthened by noting that m ^ 1 and m ^ 1. 



6.2. Collisional Dissipation of Alfven Waves 

The (field free) mean free path for Coulomb collisions in a hot accretion flow is a 
significant fraction of the local radius. This might suggest that dissipation of waves by 
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microscopic viscosity could be important. The velocity fluctuation of an Alfven wave is, 
however, perpendicular to the local magnetic field and so only cross field components of 
the viscosity tensor are important (GS). These are smaller than the field free component 
by a factor of ~ {I'pp/flpY <C 1, making microscopic viscous dissipation unimportant. 
More concretely, the dissipation range of a turbulent cascade set by microscopic viscous 
effects occurs at a scale ~ R/Re^^'^, where the Reynold's number is Re ~ vaR/ (J'± 
and /i_L ~ 0.1n/(Tp/^i?^) cm~^ s~^ is the cross field kinematic viscosity coefficient (e.g., 
Spitzer 1961). Using equation (|16D we find that Re ~ lO^V^^^l/? + 1)"^ 

mr ' and so 

rjpp ~ 10-^a-i/2/3-^/«(/3 + l f^m-'/'m'/\. (24) 

Microscopic viscous effects thus become important only on negligibly small scales. An 
analogous result holds for thermal conductivity. 

Finite electrical resistivity also leads to wave damping, which becomes important 
at a scale r^^ ~ R/ (Rem)^^^, where Rem ~ vaR/Vc is the Magnetic Reynold's number at 
the outer scale and rje ^ 10^^/Tg^/^ ~ Tg^^^ cm^ s^^ is the electrical resistivity (Spitzer 
1961). Using equation (^) yields Rem ~ 10^^(3~^^'^mr~^^'^Tg^'^ and so resistive damping 
occurs when 

r,Jpp ~ 10-^/?-l/«m-l/^ml/Vl/s^9-'/^ (25) 

CoUisionless dissipation sets in on scales ~ pp, which is ^ than the coUisional dissipation 
scales considered in this section; as the energy cascades to small wavelengths, it will 
therefore first encounter coUisionless dissipation processes. 



6.3. Implications of Particle Heating for ADAFs 

It is usual in ADAF models to specify the relative heating of protons and electrons 
by a parameter 6, the fraction of the viscous energy which heats the electrons. In order 
for the optically thin ADAF formalism to be relevant to an accretion flow, one of two 
physical situations must occur: 

1. 6 ^ 0.5 and m ^ a^. In this case a significant fraction of the viscous energy 
is transferred to the protons; by virtue of the low accretion rate (which implies low 
densities; see eq. 0]), Coulomb collisions are too inefficient to transfer this energy to 
the electrons in the inflow time of the gas (Rees et al. 1982). Consequently, a fraction 
~ 1 — 5 of the viscous energy is advected, by the protons, onto the central object. 

2. 5 ~ 1 and m ^ 10~^q;^. In this case most of the viscous energy heats the electrons; 
the flow can be advection dominated only if the electron cooling time (the time for the 
electrons to radiate their thermal energy) is longer than the inflow time of the gas. For 
synchrotron cooling, the most efficient cooling mechanism in ADAFs at low accretion 
rates (Narayan & Yi 1995b), this occurs for m ^ 10~^a'^ (Mahadevan & Quataert 1997). 
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The contribution to 6 from the coUisionless dissipation of Alfven waves considered 
in §3 is (eq. 0) 

For accretion flows with Tp ^ T^, this contribution is <^ 1. For most 50%) of the 
energy in Alfvenic turbulence to be dissipated by the mechanisms considered in this 
paper, however, the accretion flow must have (3^5 (§5.2), which is different from the 
value of /3 ~ 1 usually used in ADAF models. Furthermore, numerical simulations of 
MHD turbulence in thin disks suggest that a and P are coupled, with a ~ 0.5/(/3 + 1) 
(Hawley, Gammie, & Balbus 1996; Table 4). In this case (3^5 would correspond to 
a ^ 0.1, which is smaller than the value of a ~ 0.25 usually used in ADAF models. 
This a — (3 relationship may not, however, be applicable to ADAFs, which are radially 
convective and thus have a purely hydrodynamic source of angular momentum transport. 

ADAF models typically take 6 ~ rrie/mp ~ 10~^ (e.g., Narayan et al. 1997), but 
they are relatively insensitive to 5 so long as 5 ^ 1. The reason is that, in this limit, 
viscous heating of electrons is not their dominant heating mechanism. This can be seen 
by considering the energy equation for electrons in a hot accretion flow, 

ds de ^ , , 

a- = kTv — 
= q^^ + g'', 

where s is the entropy of the electrons per unit mass of the gas, e is the internal energy 
of the electrons per unit mass, q'^ is the compressive heating (or cooling) rate per unit 
volume, and q~ is the energy loss due to radiative cooling. The total external heating 
of the electrons, g*^"*", is a sum of the heating via Coulomb collisions with the hotter 
protons, g"^, and direct viscous heating, q'". 

For accretion rates such that 

m ^ 10-W (^^) T^/\ (28) 

Coulomb heating of the electrons dominates over viscous heating (Mahadevan 1997), 
and so the precise value of S is unimportant. The ratio of compressive to viscous heating 
is given roughly by (Mahadevan & Quataert 1997) 

gVg^ ~ lO-^r^Tgr. (29) 



For small 5, compressive heating of the electrons is more important than viscous heating, 
again making the precise value of 6 unimportant. For most of the systems to which 
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ADAF models have been applied, rh is sufficiently high and compressive heating is 
sufficiently important that only if 6 were ^ 0.03 in the inner regions of the accretion flow 
(r ~ 3 — 100, where the observed radiation originates) would the models be significantly 
modified.[^ 



6.4. Isotropy? 

The mechanisms for dissipating Alfvenic turbulence considered in this paper heat 
only the component of a particle's energy which is parallel to the local magnetic field 
(§4). In our calculations, however, we have assumed an isotropic distribution function, 
which must be justified. A standard isotropization mechanism for both electrons and 
protons is pitch angle scattering. This can occur much more efficiently than particle 
heating if the parallel phase speed of the waves responsible for the scattering is much 
less than the particle's speed (Melrose 1980); in this case isotropy could be maintained 
without significant heating. While this is a plausible mechanism for maintaining isotropy 
in hot accretion fiows, we have not investigated this mechanism in detail. 

Rather, we wish to point out a mechanism, unique to ADAFs, which maintains 
rough isotropy in the proton distribution function provided that the turbulent heating 
mechanism affects only the parallel component of the proton's energy. In addition to 
heating by the viscously generated energy, particles are heated by compression as they 
accrete inwards. Since the particles are tied to the field lines, which are being com- 
pressed, adiabatic invariance of the particle's magnetic moment requires E± oc B, where 
E± is the particle's perpendicular energy. For protons in an ADAF, the viscous heating 
rate is comparable to the compressional heating rate; this is because most of the viscous 
energy is stored as thermal energy of the protons and so all terms in the proton entropy 
equation are of the same order. If the viscous heating mechanism heats primarily the 
parallel energy of the protons, as the wave-particle interactions considered in this paper 
do, then the rough equality of viscous and compressional heating rates implies a rough 
equality of parallel and perpendicular heating, so that isotropy is maintained. Further- 
more, Mahadevan & Quataert (1997) have shown that adiabatic compression maintains 
a thermal distribution of particles even in a coUisionless gas (provided the particles are 
not compressed to relativistic energies). Therefore, given a thermal distribution of pro- 
tons at large radii in an accretion fiow, the perpendicular component of the distribution 
function remains thermal as they accrete onto the central object. 



^•^An important exception to this is the Narayan et al. (1997) model of Sagittarius A*. For this 
system, the data is so good and the estimated accretion rate is so low (m 10""*) that increasing S to 
^ 10~^ is problematic (see Fig. 4 of their paper). 
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This mechanism will not work for electrons since there is no necessary relationship 
between the compressive and viscous heating rates (the electrons are not, in general, 
advection dominated). For accretion rates satisfying equation (^), however. Coulomb 
collisions can maintain isotropy in the electron distribution function. This is because the 
characteristic timescale on which Alfvenic turbulence modifies the proton and electron 
distribution functions is ~ cu,^* .0 Since the electrons (but not the protons) are effectively 
coUisional on this timescale, isotropy is maintained. 



6.5. Comments on related work 

In a recent paper, Bisnovatyi-Kogan & Lovelace (1997; hereafter BL) suggest a 
mechanism which they feel leads to preferential electron heating in hot accretion flows. 
It is worthwhile examining this in some detail. 

Throughout this paper we have argued that, on a microscopic level, particle heating 
occurs when the MHD turbulence generated at large scales in the accretion flow cas- 
cades to small scales and is dissipated by coUisionless effects. BL seem to disregard this 
possibility, instead claiming that because the dissipation scale set by the hydrodynamic 
or magnetic Reynold's number is extremely small, plasma instabilities set in which dis- 
sipate the turbulent energy. CoUisionless effects, however, become important on length 
scales well above those set by the hydrodynamic or magnetic Reynold's number (§6.2), 
obviating the need for an appeal to plasma instabilities. 

BL argue that, on a microscopic level, particle heating in ADAFs is due to electric 
fields parallel to the local magnetic field accelerating particles to runaway velocities. By 
virtue of their smaller mass, electrons are more efficiently accelerated. In MHD, however, 
the electric field is given by E = — v x B/c + r/e V x B/ Aire. In a highly conducting 
plasma, the magnitude of a typical electric field is ~ vB/c. Furthermore, the electric 
field is primarily perpendicular to the local magnetic field and thus unimportant for 
accelerating particles. In MHD parallel electric fields arise only from finite resistivity 
corrections and are ~ vB/cRcm'-, since Rcm ^ 1 they are far too small to significantly 
accelerate particles in hot accretion flows. 

Following Bisnovatyi-Kogan & Ruzmaikin (1976), however, BL argue that, because 
the flow is turbulent, one should use a "turbulent resistivity" (of order the Shakura- 
Sunyaev turbulent viscosity) instead of the usual microscopic resistivity; in this case, 
they argue, E\\ ~ vB/c and parallel electric fields can significantly accelerate particles. 



^"^This can be estimated from the quasi-linear diffusion equations, using the results of §5.1 (i.e., it 
is only the parallel fields, which arise from kinetic corrections to GS's MHD cascade, which lead to 
diffusion of the particles in velocity space). 
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Our concern with this analysis is that, while the notion of a turbulent transport co- 
efficient may be useful in describing the large scale properties of the flow (the global 
transport of angular momentum, for example), it should not be used on a microscopic 
level. When considering microscopic processes such as particle acceleration, the global 
transport processes in the flow are unimportant. 

BL's electron heating mechanism amounts to the claim that, on a microscopic level, 
the local electric field is significantly aligned with the local magnetic field. We see no rea- 
son, however, why the MHD results described above (the local electric field perpendicular 
to the local magnetic field) should be inapplicable to hot accretion fiows. 



7. Summciry and Discussion 

Hot accretion flows, such as advection dominated accretion flows (ADAFs), are ef- 
fectively coUisionless for all but the largest scale, lowest frequency, motions of the plasma 
(§6.1). Collective plasma effects are thus hkely to be important in these systems. This 
is particularly true for particle heating since microscopic viscosity, thermal conductivity, 

and electrical resistivity arc important only at extremely small scales. Particle heating 
influences both the global structure of the flow (e.g., by determining if the protons or 
electrons are heated) as well as the observed radiation (e.g., by determining the particle 
distribution functions) . 

Particle heating is a particularly important issue for ADAF models. All ADAF 
models which have been applied to observed systems assume that the viscously generated 
energy primarily heats the protons. This enables the accretion flow to be advection 
dominated so long as the timescalc for electrons and protons to exchange energy by 
Coulomb collisions is longer than the inflow time of the gas. If viscosity only heats the 
electrons, an (optically thin) accretion flow can be advection dominated only when the 
electron cooling time is longer than the inflow time, which occurs at such low accretion 
rates that ADAF models would probably be less relevant for observed systems. Recent 
work by Bisnovatyi-Kogan & Lovelace (1997) suggests that all of the viscous energy 
in hot accretion flows does heat the electrons. In §6.5 we have argued against this 
conclusion. 

In this paper we have focused on one aspect of the particle heating problem, namely 
particle heating by the linear coUisionless dissipation of MHD, in particular Alfven, 
waves. We flnd it likely that the viscous energy in hot accretion flows resides primarily 
in MHD turbulence, making this mechanism of particular importance. An investigation 
similar to ours, but with somewhat different conclusions, was carried out independently 
by Gruzinov (1997). 
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Long wavelength Alfven waves excited on the outer scale in hot accretion flows 
are undamped by linear coUisionless effects. The nature of Alfvenic turbulence is thus 
crucial for assessing Alfven wave damping and particle heating. Recent work on Alfvenic 
turbulence by Goldreich and Sridhar (1995; GS) suggests that the turbulent energy 
cascades almost entirely perpendicular to the local magnetic field (i.e., in the inertial 
range, k± 3> k^). As a result, the cyclotron resonance, which is usually thought to be 
significant in dissipating Alfvenic turbulence, is unimportant (since the wave frequencies 
are always much less than the proton cyclotron frequency). 

We have shown, using the full kinetic theory dispersion relation for linear pertur- 
bations to a plasma, that Alfven waves with frequencies much less than the proton 
cyclotron frequency, but with perpendicular wavelengths of order the Larmor radius of 
thermal protons, are damped by transit time damping (TTD) in plasmas appropriate 
to hot accretion flows (Tp ^ and /3, the ratio of the gas pressure to the magnetic 
pressure, ^ 1). TTD is due to particles being accelerated by gradients in a wave's longi- 
tudinal magnetic field perturbation, and is the magnetic analogue of Landau damping. 
For Tp ^ Te, TTD quite generally leads to most of the dissipated wave energy heating 
the protons, since they have the larger magnetic moment and so couple better to the 
wave's magnetic field perturbation (§3.3 and §3.4). For example, if Tp ~ lOOTg, as is 
the case in the interior of ADAF models, the electrons receive only ~ 2% of the Alfven 
wave energy as it is damped in a /3 ~ 1 plasma. 

The dissipation rate of the Alfven wave increases with increasing f3 since the proton's 
magnetic moment is effectively larger and the wave-particle coupling is stronger (Figure 
2; §3.3). For ^ 5, the dissipation of Alfven waves by TTD is sufficiently strong 
to convert most (^ 50%) of the Alfvenic energy to thermal energy before the MHD 
approximations utilized in the GS analysis cease to be valid. In the GS cascade, however, 
the energy travels so quickly through the inertial range that, for /3 ~ 1, TTD is not strong 
enough to dissipate the turbulent energy. In this case, the GS analysis of the turbulent 
dynamics breaks down before most of the turbulent energy is dissipated, and it is unclear 
to us how the cascade proceeds; how the Alfvenic energy is ultimately dissipated in a 
/9 ~ 1 plasma therefore remains unresolved by our work. In §5.2 we have enumerated 
(qualitatively) what we take to be the plausible possibilities. 

So long as a significant fraction 50%) of the viscous energy heats the protons, the 
ADAF formalism is capable of describing the structure of an accretion fiow (§6.3). This 
work suggests that, for f3 ^ 5, Alfvenic turbulence in hot, two temperature, accretion 
flows preferentially heats the protons. Since we expect a signiflcant fraction of the 
viscously generated energy to reside in Alfvenic turbulence, this alone can plausibly lead 
to greater than ~ 50 % of the total viscously generated energy heating the protons. 

Wave-particle interactions are often suggested as a mechanism for forming strong 
nonthermal features in the electron and proton distribution functions. In a perpendicular 
Alfvenic cascade of the kind proposed by GS, however, the wave damping is always due 
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to particles with velocities equal to the Alfven speed (§4). In linear theory, it is thus 
impossible to accelerate particles to supra-Alfvenic velocities. For /? ~ 1 and Tp ^ 
Tg, the Alfven speed is near the thermal peak of the proton and electron distribution 
functions, so that the dissipated turbulent energy does not significantly modify the 
distribution functions from a Maxwellian. Perpendicular Alfvenic turbulence is therefore 
not a plausible mechanism for producing power law features in the proton and electron 
distribution functions or for accelerating particles to relativistic energies. This may 
exclude Alfvenic turbulence as a viable mechanism for particle acceleration in solar 
flares and accretion disk corona. 
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A. Definitions of Some Oft-used Quantities 



Quantity 


Definition 


Meaning 






mode frequency 


m 

T 


O IT) 1 \ 

2 7r/Re(u;) 


mode period 


7 


T / \ 


11 • J 

mode dampmg rate 


k 




mode wavevector 


i Ip 


qB/mpC 


proton cyclotron irequency 


Vts 




proton (s = p) and electron (s = e) thermal speeds 


Pp 


Vtp/^p 


Larmor radius of thermal protons 


Xp 


0.5klpl 


dimensionless perpendicular wavevector 


V 


kzPp 


dimensionless parallel wavevector 


a 




Shakura-Sunyaev viscosity parameter 


P 


87mkB{Tp + Te)/B'^ 


ratio of gas to magnetic pressure 


Ps 




proton (s = p) and electron (s = e) heating rates 


5 




fraction of viscous energy heating electrons 



B. The Balbus-Hawley instability in a collisionless gas 

There are three points which suggest to us that the MHD instability of Balbus & 
Hawley (BH) should apply to collisionless systems. While suggestive, these do not, of 
course, constitute a proof. The primary caveat to these comments is that, if the particle 
distribution functions are highly nonthermal, all bets are off. 

1. Perhaps the primary concern in passing from the collisional to the collisionless ver- 
sion of an instability is that collisionless dissipation mechanisms may inhibit the insta- 
bility. To linear order the axisymmetric version of the BH instability is, however, non- 
compressive and Alfvenic in character. As discussed in §3.1, in the MHD limit Alfven 
waves are undamped by linear collisionless effects, which suggests that the instability 
should not be inhibited. 

2. The collisionless limit entails the infinite conductivity limit used in ideal MHD; finite 
resistivity effects are particularly unimportant in a collisionless plasma. 

3. The Keplerian rotation frequency, which is the characteristic growth rate of the 
instability, is fio — 7 x 10^m~^r~^/^ rad s~^. The smallest characteristic frequency in a 
plasma is typically the proton cyclotron frequency, which is given in equation ([T6|) . For 
P <C 10^^, Qo ^ and so the particles are tied to the field lines, which is a requirement 
for the instability to function. In the limit of Qq ~ ^p it is unlikely that the instability 
will persist. This corresponds, however, to exceedingly small magnetic field strengths 
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Figure Captions. 

Figure 1: Properties of the Alfven wave as a function of Xp — 0.5k'j_Pp for a; <^ Q^; 
various proton to electron temperature ratios, Tp/Tg, are considered for a plasma with 
equal gas and magnetic pressure {P — 1). (a) The parallel phase speed in units of 
the Alfven speed. For Tp ^ lOTg, is nearly the same as for Tp = lO^Tg. (b) The 
dissipation per mode period, 7T. (c) The dimensionless proton (solid line) and electron 
(dotted line) heating rates (Ps)- Tp/T^ is shown along side each curve. For Tp ^ lOT,,, 
the proton heating rate, Pp, is nearly identical to the Tp = lO^Tg case, (d) The relative 
proton and electron heating rates, Pp/Pe- 

Figure 2: Properties of the Alfven wave as a function of Xp = 0.5k'j_Pp for u <^ Qp] 
various ratios of gas pressure to magnetic pressure {(3) are considered for a Tp = lOOTg 
plasma, (a) The dissipation per mode period, 7T. (b) The relative proton and electron 

heating rates, Pp/P^. 

Figure 3: The relative contribution of transit time damping {Xyy\Ey\'^) Landau 
damping {xtzl-^zl'^) to the (a) proton and (b) electron heating rates as a function of the 
proton to electron temperature ratio (Tp/Te); several f3 are considered, taking Xp = 0.1. 

Figure 4: The fractional dissipation per mode period, 7T, in the MHD limit for (a) the 
fast mode and (b) the slow mode as a function of the angle between the wavevector and 
the background magnetic field {6); various Tp/T^ are considered for a = 1 plasma. For 
the slow mode, the curves are nearly vertical displacements of each other, so fewer Tp/T^. 
are shown. 
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